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^ H 1 Abstract 

^ r*| Using first principles methods we explore the anisotropy of the spin relaxation and trans- 

verse transport properties in bulk metals with respect to the direction of the spin quantization 
axis in paramagnets or of the spontaneous magnetization in ferromagnets. Owing to the pres- 
ence of the spin-orbit interaction the orbital and spin character of the Bloch states depends 
sensitively on the orientation of the spins relative to the crystal axes. This leads to drastic 
changes in quantities which rely on interband mixing induced by the spin-orbit interaction. 
The anisotropy is particularly striking for quantities which exhibit spiky and irregular distri- 
bution in the Brillouin zone, such as the spin-mixing parameter or the Berry curvature of the 
electronic states. We demonstrate this for three cases: (i) the Elliott- Yafet spin-relaxation 

i_i mechanism in paramagnets with structural inversion symmetry; (ii) the intrinsic anoma- 

lous Hall effect in ferromagnets; and (iii) the spin Hall effect in paramagnets. We discuss 
^ the consequences of the pronounced anisotropic behavior displayed by these properties for 

spin-polarized transport applications. 

in 
q 

\Q 1 Introduction 

O 

<N 

Phenomena belonging to the field of spintronics are associated with the spin of electrons, which 
do the job of carrying information accross a device. In such a situation, the fact that the spin 
and orbital degrees of freedom of Bloch electrons in a solid are fundamentally related due to the 
presence of spin-orbit interaction (SOI) becomes of great importance. Normally, the spin-orbit 
interaction can be considered as a small perturbation compared to the other relevant energy 
scales for electrons in a crystal (such as band gaps, band widths, or exchange splittings). Its 
influence is to mix the spin and orbital character of the Bloch states for each Bloch momentum 
in a non-trivial fashion. If we consider now a non-equilibrium situation of an electron moving 
in one of the Bloch bands across the crystal under the influence of an external electric field, 
the spin-orbit mediated interaction with other Bloch states will determine its spin and orbital 
dynamics. Consider the case of a paramagnetic crystal. If we manage to make our initial 
incoming electron spin-polarized — a typical situation in a spin-injection experiment — this spin 
polarization will decrease in time due to random scattering off impurities or phonons pQ, until it 



1 



completely vanishes. Its exponential decay in time is characterized by the spin relaxation time, 
which serves as one of the most basic material parameters in spintronics [21 [3]. The presence 
of the SOI modifies in a subtle way the dynamics of Bloch electrons under an applied electric 
field, by adding a spin-dependent transverse component to the velocity. This will result in an 
anomalous Hall effect (AHE) in ferromagnets [4J, and spin Hall effect (SHE) in paramagnets [5]. 
Conceptually, understanding of these two phenomena over the past 12—15 years generated a 
cascade of novel paradigms in modern spintronics and solid state physics in general. While 
practically the AHE and SHE entered an everyday toolkit in experimental spintronics, further 
exciting research in this field is still ahead of us. In particular, a lot remains to be done concerning 
the microscopic understanding and first principles description of spin-relaxation phenomena and 
transverse transport properties in real materials. 

The crystal field in a solid is manifestly anisotropic and it lifts the degeneracy between the states 
with different magnetic quantum numbers. This results in a strong dependence of the spin and 
orbital character of the Bloch states on the choice of the spin quantization axis (SQA) or the 
direction of the magnetization in the crystal, since the matrix elements of the orbital angular mo- 
mentum operator are strongly anisotropic themselves. In ferromagnets, the crystal field splitting 
combined with the anisotropy of the orbital angular momentum operator results in a dependence 
of the eigenvalues on the direction of the magnetization, and leads to the magneto-crystalline 
anisotropy energy (MAE) [Tj, [8J — one of the most fundamental characteristics of magnetic ma- 
terials. In the field of transport phenomena in metals, the anisotropy of the electronic structure 
with respect to the magnetization direction leads to such prominent phenomena as anisotropic 
magnetoresistance (AMR) [9], tunneling anisotropic magnetoresistance (TAMR) [TOJ [H] and 
ballistic anisotropic magnetoresistance (BAMR) [12]. As in the case of the MAE, these effects 
can be already captured in many cases by considering only the changes in the band topology 
in the Brillouin zone, as the orientation of the magnetization is varied jT3l [T2j HH [TO]. This 
situation is in contrast to the case of the AHE and SHE, which are often governed by band 
degeneracies at the Fermi level [15J. In this case the dependence of the eigenenergies on the 
SQ A/magnetization is either absent or can often be neglected, while the anisotropy of the spin 
and orbital resolution of the wavefunctions becomes of primary importance, and could lead to 
very large values of the anisotropy of the Hall conductivities, as speculated already by Fivaz in 
1969 for the anomalous Hall effect [16J. 

The significant anisotropy of the spin-relaxation and Hall effects is a valuable tool for tuning the 
transport properties of spintronics devices. Since such anistropy is an intrinsic property of the 
mono-crystalline solid, it should be properly averaged when using polycrystalline samples, as well 
as when considering the effect of temperature and magnetization dynamics on the measured spin- 
polarization or transverse current |17j . Experimentally, only the anisotropy of the anomalous 
Hall effect has been researched in the past and in many cases a very large anisotropy was 
found [IHl [IB [201 Ell [221 ESI [2^ while evidence of anisotropy in the spin Hall effect [26] and 
spin-relaxation has been presented only recently [27] [28] . 

Here, we review the present theoretical understanding of anisotropy in three phenomena occuring 
in perfect crystals: (i) spin relaxation, (ii) intrinsic anomalous Hall effect, and (iii) intrinsic spin 
Hall effect, Fig. 1. We will focus on the developments which took place over the past years jTTl 
[29l [301 EU |32l [331 El]- In particular, we present agruments and show from first-principles 
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calculations that due to the sensitivity of spin-relaxation and Hall effects to the SOI-mediated 
coupling between (nearly) degenerate states in the vicinity of the Fermi level, the anisotropy 
of these effects can be gigantic, and has in principal no theoretical limit. Manifestly, for some 
directions of the SQA and magnetization in the crystal the spin relaxation rates and Hall currents 
can be suppressed by orders of magnitude, or even display a change of sign in corresponding 
conductivity components. We discuss possible applications of such large anisotropics, encourage 
further experimental studies in this area, and emphasize that a wide range of materials exhibit 
anisotropic transverse transport and spin relaxation, from bulk solids to surfaces and interfaces 
with essentially lowered lattice symmetry. 

1.1 Spin relaxation in paramagnets 

As a first example of a situation in which the importance of the anisotropy of the wavefunctions 
with respect to the choice of the spin quantization axis is very pronounced we consider the 
phenomenon of spin relaxation. To be concrete, here we concentrate on the Elliott- Yafet spin- 
relaxation mechanism, dominant in materials with structural bulk inversion symmetry [TJ [3], 
which is due to scattering of electrons off phonons or impurities. Owing to the presence of 
spin-orbit coupling (SOC) in the system such scattering events will flip the spin of the electron 
with a certain probability, which depends on both the wavefunctions of the ideal crystal and 
the scattering potential. However, according to the Elliott approximation p], an estimate of 
the ratio between momentum- and spin-relaxation times, r p and Ti, can be given in a first 
approximation by neglecting the form of the scattering potential as follows: t v /T\ « 46 2 , where 
b 2 is the Elliott- Yafet parameter (EYP) defined below, which is an intrinsic property of the ideal 
crystal [35J. 

The coexistence of time-reversal and space-inversion symmetries implies that the eigenenergies 
of the crystal at any given Bloch momentum k are at least two-fold degenerate. Following 
Elliott, we write the corresponding states as 

Vt(r) = [«ks(r) |t)s + Mr) IDs] e ikr , (1) 
V4(r) = [«- ks (r) |l)s-&_ ks (r) |t)s] e lk ' r . (2) 

The two spin states |t)s an d \\)s are eigenstates of s-S, where s is the unit vector along the chosen 
SQA, S = |cr^ is the spin angular momentum operator, and cr p are the Pauli matrices. So, for 
example, \^) z and \l) z are the eigenstates of the S z operator. More generally, the reference frame 
is specified by the SQA direction s, which is chosen to coincide with the polarization direction 
of the initial/injected spin population. The functions aks(r) and 6ks( r ) exhibit the periodicity 
of the crystal lattice, and we define 6^ as the unit cell integral / d 3 r |6ks( r )| 2 (similarly for 
a ks> wmcn satisfies aj^ = 1 — 6^ ). 

For fixed s, the degenerate and states (and the corresponding aks(r) and &ks( r )) can 
be chosen, via a linear combination, such that the spin-expectation value Sks = I'S'sIr ) 
is maximal. The spin mixing parameter is then given by 6^ = 1/2 — Skg/h, and is usually 
small, due to the weakness of the SOC. In this case the Bloch states are of nearly pure spin 

character. (Thus, aks represents the "large" component of the spinor, while &ks the "small" 

"t i 

component. The relation between the large or small components of tp^ and -0^, = &lks 
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and &ks — ^_ks' f°H° ws from space and time inversion symmetry). However, at special spin- flip 
hot-spot points in the band structure, e.g. accidental degeneracies, Brillouin zone boundaries or 
other high symmetry points [35] . 6^ may increase significantly up to |, which corresponds to 
the case of fully spin-mixed states. Generally, the distribution of the spin-mixing parameter for 
a metal with a complicated Fermi surface can be very non-trivial. The Fermi-surface averaged 
Elliott- Yafet parameter is given by 

6? - 1 - [ |6kg|2 d 2 k (3) 

where v j p(k) is the Fermi velocity. The normalization by the density of states at the Fermi level, 
n(Ejr) = 1/h J FS ^^(k)! -1 d 2 k, ensures that < bj < \. 

For the ensuing discussion it will be useful to divide the spin-orbit operator into spin-conserving 
and spin- flip parts, £(LS^) and £(LS^), given respectively by the first and second terms of the 
following expression [30J: 

£L • S = ZL,Ss + £ (L+S7 + L7S+) /2 . (4) 

Here £ is the spin-orbit coupling strength, L is the operator of the orbital angular momentum, 
Lg = L • s, Sq = S • s, and and are the raising and lowering operators for orbital and spin 
angular momenta. Acting on a state of the crystal obtained without SOC, the spin- flip part of 
the SOI can flip its spin, while the spin-conserving part will keep it intact. It is clear that the 
dot product L • S is independent of s, leaving the eigenenergies of the Hamiltonian invariant. 
However, the spin-conserving and spin-flip parts, separately, depend on the choice of the SQA. 
In ferromagnets, the spin-conserving part of SOI is the one which is largely responsible for the 
values of the magneto-crystalline anisotropy energy and orbital moments [HI E] • In paramagnets, 
the Elliott-Yafet spin-relaxation mechanism is driven by the spin-flip part of the SOI. 

In an experiment, the spin polarization of the electrons subject to spin relaxation is defined by 
the direction of the external magnetic field (e.g. in conduction-electron spin resonance experi- 
ments) or by the polarization of ferromagnetic leads (e.g. in spin-injection experiments). In a 
paramagnet, the choice of the spin-quantization axis, determined by the direction of the spin 
polarization, does not influence the band energies, and its most important manifestation is in the 
changes of the orbital and spin character of the Bloch states. Experimentally, the dependence 
of the spin-mixing parameter on the SQA was observed in supported graphene layers [27] and 
in semiconductors [28]. However, no microscopic theory of anisotropic spin-relaxation which 
explicitly refers to the anisotropy of the Bloch states has been given. In bulk metals with inver- 
sion symmetry, which are at the focus of this article, the Elliott-Yafet mechanism is dominant. 
Below we will demonstrate that indeed the anisotropy of the EYP in metals can be gigantic. 

1.2 Anomalous Hall Effect in ferromagnets 

A second phenomenon, for which multiple experimental studies of anisotropy exist [HI \19\ [20J, 
EU [22j [231 EH EH], but no quantitative theoretical argumentation for its emergence had been 
presented until recently [T7j , is the anomalous Hall effect [36J, [371 EE]- The essence of the 
AHE in a ferromagnet lies in the generation of a charge current J# transverse to the electric 
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Figure 1: Setup of the anisotropic spin-relaxation and transverse spin and anomalous Hall 
effects in metals. Anisotropic spin-relaxation (left): an electron with a certain direction of spin, 
injected into a solid which exhibits anisotropic spin-relaxation will lose the memory of its initial 
spin polarization over a period of time, which depends on the direction of the spin-polarization 
in real space (z or x, red and blue arrows correspondingly). Anisotropic spin Hal effect (middle): 
the magnitude of the spin current and the direction of its spin polarization (related to the sign 
and magnitude of the spin accumulation on the surfaces of the sample) measured along a certain 
direction in a crystal depends on the direction of applied electric field E which generates the 
longitudinal electric current Jo- Anisotropic anomalous Hall effect (right): for a fixed direction 
of applied electric field E and corresponding longitudinal electric current Jo (out of the plane of 
the manuscript), the direction and magnitude of the Hall current J# depend on the direction 
of magnetization M in the sample. Note that for a general direction of M away from high- 
symmetry axes in the crystal, J# can be non-orthogonal to M. 



field E (and corresponding "diagonal" current Jo), without any applied magnetic field fT5] . 
Phenomenologically, the relation between the ith component of the Hall current and the jth 
component of the electric field is the following: 

Jh,i = VijEj, (5) 

where &ij are the components of the anomalous Hall conductivity (AHC) tensor. Since in a 
ferromagnet the anomalous Hall conductivity tensor is second rank antisymmetric, the AHC 
tensor can be also seen as the AHC vector cr, whose components are related to the components 
of the AHC tensor as G{ — \ Y^jk tijk&jk, through the Levi-Civita tensor e^: 

J H = Exa. (6) 

In general, there can be different contributions to the AHC in a ferromagnet. In a real material 
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contributions to the anomalous Hall current which originate from scattering of electrons off 
impurities or due to disorder are always present — this is the so-called extrinsic AHE [39 , 40 , 33J. 
The second important part of the AHE signal — which comes solely from the electronic structure 
of the pristine crystal — is the so-called intrinsic AHE [JT]. Namely, the SOC in a perfect crystal 
gives rise to a transverse spin-dependent "anomalous velocity" of electrons propagating along the 
direction of the applied electric field — thus leading to the intrinsic anomalous Hall current [UH2]. 
In the present work we will focus exclusively on the intrinsic contribution to the AHE. 

The intrinsic anomalous Hall conductivity is determined by the electronic structure of the pris- 
tine crystal, which can be accurately calculated using modern first principles methods, see for 
example a recent review by Gradhand et al. |13] . Several investigations of intrinsic AHC using 
the first principles methods have been done, for instance, in SrRuOs [44l[45], Fe [21 [46], Ni [47] . 
Co [HI [T7] and other ferromagnets. For these materials, the calculated intrinsic AHC agrees 
well with the experimental values, except for the case of fee Ni [47], which is most probably 
due to effect of electronic correlations [33J HE]. It is therefore a common belief that the AHE 
in moderately resistive samples of itinerant ferromagnets is often dominated by the intrinsic 
contribution. 

The intrinsic AHC considered in this work can be obtained via the linear response Kubo formula 
for the off-diagonal components of the conductivity tensor a: 



which relates the conductivity tensor to the Brillouin zone (BZ) integral of the fc-dependent 
Berry curvature tensor Q. In the latter expression and are respectively the occupied 
(o) and empty (e) one-electron spinor Bloch eigenstates of the crystal, e n \^ and £ m k are their 
eigenenergies, and v\ and vj are the Cartesian components of the velocity operator v. The 
Berry curvature ft appearing in the equation above is the very same quantity which arises when 
the adiabatic dynamics of electrons in the reciprocal space is considered [12]. In particular, 
the Berry phase acquired by a Bloch electron as it traverses a closed path in the Brillouin 
zone can be calculated as an integral of the Berry curvature over the enclosed area [42[ [49] . 
Mathematically, the Berry curvature is the curvature of the fibre bundle of occupied electronic 
states in an insulator, and its integral over the whole torus of allowed Bloch vectors provides the 
value of the quantized transverse charge conductivity, as first demonstrated for the case of the 
quantum Hall effect by Thouless et al. [50] . The Berry curvature is a key quantity in the field 
of Chern and topological insulators [HI [52]. Overall, Eq. Q manifests the topological nature 
of the intrinsic anomalous Hall effect in metals. 

Let us briefly outline the concept of anisotropy as it applies to the AHE. The definition is slightly 
more complicated than for the Elliott- Yafet parameter discussed previously, owing to the vector 
nature of the anomalous Hall conductivity. The parameter with respect to which the anisotropy 
of the AHE is studied is the direction of the magnetization M in the crystal. The anisotropy 
of the AHC with respect to M is two-fold: not only the magnitude of a depends on M, but 
also the direction of cr displays a non-trivial dependence on the direction of magnetization. For 
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a high symmetry direction of M the AHC vector is aligned with the magnetization so that the 
Hall current is perpendicular to it. For a general direction of M away from high symmetry axes 
in the crystal the AHC vector can deviate from the direction of M, in which case [20| [T7] : 



where <T||(M) is aligned with M while a±(M.) is perpendicular to it [T7] (see also Fig. 1). 
The microscopic origin of the AHE anisotropy is clear from the expression Q for the Berry 
curvature, according to which both dependence of wavefunctions as well as eigenenergies on the 
magnetization direction leads to the anisotropy of the AHC. It is important to realize that in 
contrast to the case of paramagnets with inversion symmetry (considered in the following with 
respect to the anisotropy of the SHE and EYP), for which also the eigenspectrum does not 
change with the SQA, the dependence of the wavefunctions on the magnetization direction in a 
ferromagnet is far more complex, owing to broken time-reversal symmetry. Also, the anisotropy 
of the velocity matrix elements has to be taken into account in uniaxial crystals. The largest 
part of this work is dedicated to analyzing the anisotropy of the AHE in uniaxial crystals. 

1.3 Spin Hall Effect in paramagnets 

The spin Hall effect in paramagnets consists in generation of a spin current orthogonal to the 
direction of an applied electric field E. In a simple picture, the spin current in the SHE can be 
seen as two anomalous Hall currents, propagating in opposite directions for spin-up and spin- 
down electrons. In contrast to the AHE, where the direction of the Hall current is uniquely 
determined by the directions of E and M, the spin Hall current propagates in all directions 
othogonal to E. For each of the directions of the spin current, the "physical" spin-quantization 
axis is determined by the direction of the current's spin polarization. First proposed theoretically 
in 1971 [5], the SHE was "re-discovered" in 1999 [53] . and eventually experimentally observed in 
2004 [53], triggering development of new directions in spintronics [55 [ 56J and further research in 
the direction of quantum spin Hall insulators [571 EH1 EH1 EQ] . In analogy to the anomalous Hall 
effect, the observed SHE in metals contains two types of contributions: one extrinsic (driven by 
disorder), and the other intrinsic (diorder-independent) [391 ED E3 E2] • And while very often 
the spin Hall effect is associated with the resulting spin accumulation at the boundaries of the 
sample, employing inverse SHE it is possible to measure directly the spin Hall conductivities 
(SHCs), which are much easier to treat theoretically with ah initio methods. As in the case of 
the AHE, for transition metals the experimental SHC values agree very often with the values 
obtained from first principles calculations for the intrinsic SHE. 

In the first principles calculations presented below we consider only the intrinsic [571 ESI EH EHl 
|6T] contribution to the SHC, which results from the virtual interband transitions in the presence 
of an external electric field. It may be expressed using a linear response Kubo formula analogous 



er(M) = <7ii(M) +<tj_(M), 



(8) 



to Eq. Q for the AHC: 




o,e 



(VVik | Qj | ^mk) (^mk | Vj | VVik) 



(9) 



^•(k) =-2Im^ 



n.m 
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where Qf is the spatial i- and spin s-component of the spin velocity operator, and tensor ft 
is sometimes referred to as the "spin Berry curvature". If only the spin-conserving part of 
the SOI is taken into account, the spin-projection along the direction of the spin-polarization 
of the current s is a good quantum number, and the spin velocity operator may be written as 
Ql = | Vi(j v s . In this case the SHC equals twice the value of (scaled) anomalous Hall conductivity 
for spin-up electrons only. Here, a v s is a Pauli matrix used to express the s-component of the 
spin operator. In order to treat the spin-nonconserving part of the SOI correctly, we used the 
definition of the spin current density operator given in Ref. [65 j . 

To our knowledge, experimentally, the anisotropy of the spin Hall effect has been discussed only 
once for AlGaAs quantum wells [26]. In metals the anisotropy of the SHE was investigated 
recently by Freimuth et al from first principles [29] . In many aspects, the SHC anisotropy is 
analogous to that of the anomalous Hall conductivity. It is remarkable, however, that due to 
the higher symmetry of the problem the anisotropy of the SHC in transition metals is what we 
call purely geometrical. By this term we mean that it is exactly absent in case of a cubic crystal, 
while generally the dependence of the magnitude and spin polarization of the spin current on its 
direction can be reconstructed exactly from corresponding values for high-symmetry directions 
in the crystal. This is in sharp contrast to the behavior of the AHC or EYP, which exhibit 
anisotropy already in cubic crystals, and for which the dependence of the magnitude of the EYP 
(anomalous Hall current) on the direction of the SQA (magnetization) cannot be reconstructed 
from the respective "high-symmetry" values. The middle part of this review is dedicated to the 
anisotropy of spin Hall effect in transition metals. 

2 Computational methods 

For calculations of the Elliott-Yafet parameter and corresponding Fermi surfaces we used density 
functional theory in the local density approximation |66j to calculate the underlying electronic 
structure. For the self-consistent calculations we employed the Korringa-Kohn-Rostoker (KKR) 
Green-function method [67J in the atomic sphere approximation and solve the Dirac equation 
with angular-momentum expansion up to ^ max = 4. The Fermi surface is determined by the 
KKR secular equation, det(./W(k, Ep)) = 0, which is equivalent to the condition that at least one 
eigenvalue of the KKR matrix A4 vanishes. We search for the k-vectors fulfilling this condition 
with a tetrahedron method using linear interpolation of the complex eigenvalues of A4. We 
choose a grid of 200 k-points for each direction in the full Brillouin zone, resulting in about 10 7 
Fermi-surface points. We followed the procedure described in Ref. |68j to maximize the spin 
component Sks at the Fermi-surface points. The integration (Eq. ([5])) is done by evaluating the 
integrand at the Fermi-surface points and interpolating linearly within the connecting triangles 
(the details of the method will be published elsewhere). 

For calculations of the intrinsic anomalous Hall and spin Hall conductivities we employed the 
full-potential linearized augmented plane- wave (FLAPW) method, as implemented in the Julich 
code FLEUR [69] . We used the generalized gradient approximation to the DFT and experimental 
lattice constants of the transition-metals. The self-consistent calculations with SOC were done 
in second variation with k max between 3.7 and 4.0 a.u. -1 and about 8000—16000 fc-points in 
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the full Brillouin zone. For intermediate alloys, for instance, (Feo.5Coo.5)Pt, the virtual crystal 
approximation (VCA) was applied on the 3d atomic sites, where the composition-averaged core 
potential was used instead of that of pure 3d elements, together with corresponding number 
of valence electrons, and interpolated lattice constants from the neighboring compounds [32]. 
For the calculations of the conductivities we applied the Wannier interpolation technique of 
Wang et al. [46J. We followed the method introduced in Refs. [70] and [71] to construct the 
maximally-localized Wannier functions (MLWFs) from the FLAPW Bloch states ^km : 

W Rn (r) = -j- £ e ~ lk R E U ^m(r), (10) 



k 



where Wn n denotes the n-th WF centered at lattice site R, Umn refers to the unitary transforma- 
tion among the Bloch states at k which minimizes the spread of the Wannier functions. Using 
the self-consistent charge density with SOC included, 18 spinor MLWFs per transition-metal 
atom, corresponding to d-type orbitals, were generated using wannier90 code [72] . 

Working in the basis of the maximally-localized WFs allows us to construct a real-space tight- 
binding Hamiltonian of the crystal, which can reproduce the electronic bands with any given 
accuracy at any fc-point in the Brillouin zone, given that a necessary number of fc-points was 
used for the generation of the WFs [46J. The real-space tight-binding hopping parameters can 
be calculated as: 

ffmm'(R) = ^E € kne" lk - R (u^U^,, (11) 
kn 

where i7 mm /(R) denotes the hopping parameter between Wannier orbitals Wiw(r) and Wo m (r). 
Based on those parameters, the Hamiltonian i?(k), matrix elements of the velocity operator as 
well as charge and spin Berry curvature in reciprocal space can be efficiently evaluated using 
the Wannier interpolation technique [46] . 

In the section on the anomalous Hall effect, we evaluate the derived perturbation theory ex- 
pressions for the AHC in L1q FePt. In order to apply the perturbation theory in the basis of 
Wannier functions, we use the basis of Wannier functions constructed without SOI to calculate 
the matrix elements of the spin-orbit interaction. To do this, the scalar-relativistic Hamiltonian 
without SOI is set up for the majority and minority states, and diagonalized in order to obtain 
the Bloch functions ^^(r), w ^ n a — t or I- The matrix elements of SOC in the basis of Bloch 
states can be then calculated: 



r dr 



C' > > (12) 



where L M is the atomic orbital momentum operator associated with atom \i (with the potential 
V^). In the scalar-relativistic approximation, the Hamiltonian iJ^ m ,(R) can be obtained as: 

Km' (H) = ^ E ^e- lkR (if™,) ' U^.. (13) 
kn 

Likewise, the matrix elements VL^ , are transformed into the basis set of Wannier functions: 

K^(R) = ^ E V Z^ R (uSLY^- (14) 

krm' 
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The complete Hamiltonian with SOC in the WF-basis is then given by 



(15) 



Such a separation enables us to perform the perturbation treatment of SOC. By calculating 
from first principles an atomic shell averaged SOC parameter £ d = ^ ~ dV Jr^ ) •> ^ * s possible to 
write the SOC operator approximately as £L • S, where L (S) is the total orbital (spin) angular 
momentum operator. Calculated in such a way the SOC strength £ for Pt (about 0.6 eV) is one 
order of magnitude larger than that of 3d elements, for instance, £ = 0.06 eV for Fe. 

3 Anisotropy of spin relaxation in metals 

First of all, let us work out the perturbation theory expression for the spin-mixing parameter of 
a certain state ^ n (we omit the explicit fc-dependence for the moment). Let us assume that the 
spin-conserving part of the spin-orbit interaction has been included in the Hamiltonian, which 
has as an eigenstate. In this case, since the spin-conserving SOC keeps the spin a good 
quantum number, in a paramagnet with structural inversion symmetry, the state i\) n can be 

characterized by a certain value of spin, say, ip n = ipX- This has an exact replica but of opposite 

i "t 
spin ipn. Upon including into the picture the spin- flip SOI, will acquire an admixture of the 

down spin which we will denote as (^n)^? and which corresponds to the part that includes b\& 

in Eq. It is clear that since spin- flip SOI is off-diagonal in orbital character, the down-spin 

"t i 
admixture of ipli does not come from interaction with ipn at the same energy, but comes from the 

interaction with the other states in the system. Within first order non-degenerate perturbation 

theory (^n)^ can be calculated as: 



Since the spin-mixing parameter b 2 n is equal to |(^)^| 2 , we readily obtain from the latter ex- 
pression that 



Therefore, the spin-mixing parameter of a certain state is just a sum of amplitudes for SOC- 
mediated spin-flip transitions from this state to other states and back. This picture of the 
Elliott- Yafet parameter in solids had been suggested by Elliott already in 1954 [lj. Later on in 
this work, we will apply a similar approach in order to perform a perturbation theory analysis 
of the Hall effects. 

Before proceeding with ab initio calculations, we consider a simple model which is able to 
capture the origin and essential properties of the anisotropy of the spin-mixing parameter in a 
solid. Namely, let us consider six p-orbitals, p a x , p G y and p a z with a = (t>i) standing for the spin 
of the orbit als. In order to consider the spin- mixing separately, we explicitly separate the two 
SOC terms. We have: 



H = H + Z(LS)^ + £(LS) n = diag(£, e + 5, e + A) ® l 2x2 + £(LS) n + f (18) 




(16) 




(17) 
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where in the on-site part the four states p° and p G y are chosen to be almost degenerate at energy 
£ (separated by energy 5: 5/ A « 1), and the p a z orbitals are shifted to higher energy e + A in 
order to mimic the uniaxiality of the lattice. The SOC strength is given by £, with £/A << 1. 
The energetic levels and their orbital character without spin-orbit are shown in the left column 
of Fig. [2j in which 5 was put to zero. 

Let us first consider the case when the SQA points along the z-axis and 5 = 0. When only 
spin-conserving SOC is added to i?o, the eigenstates are [yf x ± ipy)/y/2 and p a z (a =1\>0? with 
energies as sketched in Fig. [2] The only non-vanishing matrix elements of the spin-flip SOC 
are (pi — ipl\LS^\pi) = 2 and (pi + ipy\LS^\pl) — —2 which come from the states that are 
well-separated in energy. According to Eqs. ( fl6| ) and (17) this leads to a small admixture of the 
p^-state of opposite spin in the lowest lying eigenstate, and corresponding spin-mixing parameter 
of the order of (£/A) 2 when the spin- flip SOC is included. On the other hand, when the SQA is 
chosen along the x-axis, the spin-conserving part of SOC mixes small £ / A-portions of p G y with 
p a z orbitals, see Fig. [2j Now, there are four non-vanishing matrix elements of spin-flip SOI, all 
of order 1, among which two transitions are very close in energy with a separation of ~ £ 2 /A> 
giving rise to a much larger spin-mixing of the order of A/£. This results in a very strong 
spin-mixing between the two low-lying orbitals when the spin-flip SOC is added, and leads to a 
very large spin-mixing parameter of the lowest-lying state of the order of (A/£) 2 , meaning that 
higher-order perturbation theory is needed since the spin-mixing parameter cannot exceed |. 
The resulting orbital and spin character of the states when the Hamiltonian with complete SOC 
is diagonalized, is presented in the right column of Fig. [2j Note, that the final eigenenergies are 
the same, while the character of the states is different among the two directions of the SQA. 
Obviously, the resulting anisotropy of the spin-mixing parameter of the lowest-lying state with 
respect to the choice of the SQA is very large if £/A << 1. 

The model presented above allows us to make some statements concerning the general conditions 
under which a large anisotropy of the spin-mixing parameter in a metal can be expected. First of 
all, crucial is the presence of a degeneracy or near-degeneracy at Ep, of Bloch states originating 
from the atomic orbitals </> m and m /, with the orbital characters \m — m'\ ^ 1, which are the 
eigenstates of the Lg operator for some direction of the SQA (say, z). In this case {\m — m'\ ^ 1) 
no direct coupling is allowed between them by the spin-flip part of the SOC Hamiltonian, and the 
system is "protected" against large- amplitude spin- flip transitions, since the spin-mixing occurs 
due to interaction with other, energetically different, states. Correspondingly, the further away 
these other states are from the Fermi energy, the smaller the spin-mixing parameter will be, as 
exemplified for our model in Fig. [3] where at S = (we remind that S stands for the splitting 
between p x and p y orbitals), the parameter A is varied. In this case also relative position of 
the states (j) m and (f) m t with respect to each other is not that important for the spin-mixing 
parameter, see dotted line in Fig. [3j in which at constant A, parameter 8 is varied. On the 
other hand, for the SQA along the x-axis, the spin-mixing between the nearly degenerate states 
(j) m and <j) m / is favored and reaches very large values, decaying as a function of the separation 5, 
see full line in Fig. [3j The spin-mixing of the nearly degenerate states with the other states is 
minimal on the other hand, but it reduces the overall value of the spin-mixing parameter. The 
suppression of the spin-mixing parameter of the states at the Fermi energy due to interaction 
with the higher lying states is reduced, the further the latter are from the Fermi energy, as 
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Figure 2: Eigenvalues of the spin-degenerate p-states described by the Hamiltonian diag(s, £, £ + 
A) ® 1.2x2 (left column), including the spin-conserving part of SOC for two different quantization 
axes (middle column), and including both spin-conserving and spin- flip parts of SOC for two 
different quantization axes (right column). For each doubly degenerate state the orbital and 
spin decomposition is given with a label. Note that while the energy spectrum is identical for 
both SQA's in the right column, the spin character of the two low-lying states is different, which 
gives rise to the anisotropy of the spin-mixing parameter. In the labels, "e" and 'V denote a 
small admixture of the corresponding state of the order of £/A. 

clearly visible in Fig. [3j in which A is varied at constant 5 = 0. Overall, by examining Fig. |3j 
we conclude that the largest anisotropy of the spin-mixing parameter in a metal will be favored 
when the states with \m — m f \ ^ 1 at the Fermi are perfectly degenerate, and are positioned 
very far away from other states. 

Let us now turn to hep osmium, which we choose as an example of a typical transition metal 
with a uniaxial crystal structure and the properties of which we calcualted from first principles. 
First, we take a look at the bandstructure of Os along the T — A path from the Brillouin zone 
center along the z-axis, presented in Fig. [3] (right). The splitting of the two bands (each band is 
doubly degenerate) which cross the Fermi level here ("1" and "2", full lines), is due to the spin- 
orbit interaction, as can be verified from the fact that they fall on top of each other when scaling 
down the spin-orbit coupling strength. Also, without SOC, it is straightforward to determine 
the orbital character of the bands: in this case the bands have a d + i and d-i character, which 
are superimposed to form the d xz and d yz states. Overall, we have all prerequisites for a large 
anisotropy of the Elliott- Yafet parameter at this point of the Fermi surface, according to the 
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Figure 3: Left: Dependence of the spin-mixing parameter of the lowest lying state in the model 
from Fig. [2] on the separation S (in eV) between nearly degenerate p x and p y states at constant 
separation with the p z orbital A = 1 eV (left), and on the separation of degenerate p x and p y 
states (5 = eV) with the p z orbital A (in eV, right). The SOC strength of 0.1 eV was taken 
for these calculations. Clearly, the largest anisotropy of the spin-mixing parameter is acquired 
when the states p x and p y are perfectly degenerate and are lying far away from the p z orbital. 
Right: Band structure of hep osmium around the Fermi level, in the direction r — A of the BZ 
with applied B-fleld of 40 meV. The exchange splitting of the two bands crossing the Fermi 
energy (1 and 2) depends on the direction of B, reflecting the anisotropy of the spin-mixing at 
the hot-spot "H" in Fig. [4} 

arguments presented above, since bands "1" and "2" are well-separated from other bands at the 
crossing with the Fermi level. Before proceeding with an explicit calculation of the Elliott- Yafet 
parameter distribution over the whole Fermi surface of Os, we perform a numerical experiment 
in order to examine the anisotropy of the response of the bands in Fig. [3] to a small Zeeman-like 
field B with the magntiude of 40 meV. The small Zeeman field which we apply by hand lifts the 
remanent degeneracy owing to the coupling to the Bloch states of the form B • cr p , which breaks 
the time-reversal symmetry and defines a spin-quantization axis in the direction of B. In Fig. [3] 
we clearly observe a splitting of bands "1" and "2" for B along the c-axis in the crystal (dashed 
lines). However, for B in the afr-plane, the degenerate pairs "1" and "2" do not split (solid lines), 
which marks a very anisotropic response to a Zeeman magnetic field. We can understand this 
result by employing the perturbation theory arguments from before: in first order, the energy 
shift of a state due to the presence of a small Zeeman field is proportional to its spin polarization, 
which is similar to that depicted in Fig. [2| Correspondingly, while in the latter case the states 
are fully spin mixed and the bands do not split, in the former situation a large Zeeman splitting 
is achieved. 

The calculated Fermi surface of Os presented in Fig. Qa-c) consists of two nested sheets, a 
surrounding surface crossing the Brillouin zone boundary and little hole pockets (denoted by 
"P"). The latter ones are ellipsoids in an extended zone scheme centered around a point on the 
Brillouin zone boundary. Analyzing the distribution of the spin mixing parameter on the 



13 



0.1 0.2 0.3 0.4 0.5 -0.5 -0.2 0.2 0.5 




Figure 4: Fermi surfaces of Os (a-c) and W (d-e). For an illustration of the nested sheets only 
half of the Fermi surface of Os is shown. The Elliott- Yafet parameter 6^ is shown with a color 
code on the Fermi surface with the SQA along the c-axis (a) and in the a6-plane (b). Red arrows 
at the left-lower corner of (a) and (b) indicate the direction of the SQA. The difference of 6^ 
between the two directions of s is shown in (c). Analogously, for s along [001] and [111] 
in W is shown in (d) and (e), respectively. The averaged values of 6? over all directions of s, 
corresponding to poly crystalline samples, are 0.0666 for Os and 0.0627 for W. Note that the 
color scale on the right refers only to (c). Taken from [34] . 



Fermi surface, we observe a strong dependence on the SQA, evident from comparing Figs. |4](a) 
and (b). For s along the c-axis of the crystal, Fig. Q a), the spin mixing is relatively uniform 
(^ks ~ 0.05) for large areas of the Fermi surface, reaching higher values near the pockets. 
However, this picture changes drastically when s is in the afr-plane (Fig. |4^b)). In this case, the 
areas with full spin mixing (red, 6^ ^0.5) are prominent, most clearly at the caps of the two 
nested Fermi-surface sheets, indicated by "H", which are formed by bands "1" and "2" crossing 
the Fermi level in Fig. |3| Additionally, large areas with smaller, but still strong spin mixing 
(^ks ~ 0-3) are visible, e.g. in the region denoted by "B". Overall, for the two considered cases 
there is a strong qualitative difference in the k-dependent spin- mixing parameter 

As for the Fermi-surface averaged values &?, we find 6? of 4.85 x 10 -2 and 7.69 x 10 -2 for s along 
the c-axis and the a6-plane, respectively, yielding thus a gigantic anisotropy of the Elliott- Yafet 
parameter, defined as A = [maxg(6?) — min§(&|)] /ming(6?), of 59%. On the other hand, the 
anisotropy A with respect to rotations of the SQA within the a6-plane is negligible. These two 
limiting cases are depicted in Fig.JHJ^b), in which the value of b\ is shown as a function of s for all 
possible directions of s. The absent (or very small) anisotropy in the a6-plane is reflected in the 
rotationally invariant color-scale around the c-axis, as opposed to the large difference between 
the afr-plane and the c-axis. The difference of for the two limiting cases of SQA for each 
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point at the Fermi surface is shown in Fig. Qc). Large areas of the Fermi surface show small 
orientational dependence of b^. (white areas). The anisotropy at the hot spots is very large, 
but the sign is different between the hot spot "H" and the hot spots near the pockets. The 
magnitude of the effect is strongly enhanced by the large extension of the two near-degenerate, 
parallel sheets of the Fermi surface, resulting in a spin- flip "hot area" around "H" instead of a 
single "hot spot". In addition, the reduced symmetry helps: if the crystal had cubic symmetry, 
then upon change of the SQA from z to x the effects at rotationally equivalent parts of the 
Fermi surface would mutually cancel. 

Next, we analyze the hot-spot contribution to the averaged 6? and the anisotropy A. We perform 
integrals similar to Eq. Q, but restricting the integration to the part of the Fermi surface where 
b^~ lies in certain intervals, X{ < bj^ < a^+i, with X{ = 0,0.05,0.10, ... This integration results 
in values 6? which form the histogram presented in Fig. [5|a). As we can see, for the SQA 
along the c-axis, 6? is mainly determined by regions with relatively low spin-mixing parameter 
(^ks < O-l^), leading to the total value of 4.85 x 10 -2 (denoted by the solid arrow). For s in 
the afr-plane there is also a considerable contribution from regions with 6 2 g > 0.15 increasing 
the total value to 7.69 x 10 -2 (dashed arrow). Comparing the two histograms for different 
SQA, we can draw conclusions about the respective contribution of each region to the total 
anisotropy, A = (b 2 ab — b 2 )/b 2 . Interestingly, the anisotropy originates not only from the hot 
spots with > 0.35 leading to A — 12%, but mainly from the areas with smaller spin mixing 
0.15 < — 0-35 around the hot spots and regions "B", resulting in A = 49%. The larger area 
with low spin-mixing, 6^ < 0.15, does not contribute to the anisotropy significantly (A = —2%). 

Let us now turn to tungsten, which has a bcc lattice structure, and see whether higher symmetry 
of the lattice brings qualitative changes in the anisotropy of the EYP. When s || [001], 6^ exhibits 
hot-spots in directions perpendicular to s (denoted by "C") (Fig. Qd)), following a formation 
scenario similar to that at the "H" -point in Os. Additionally, many states with smaller spin 
mixing (0.2 < 6^ < 0.3) are present at the Fermi surface, leading to 6? = 6.49 x 10 -2 . For the 
SQA along another high symmetry direction of the lattice, s || [111] in Fig.Qe), the intensity at 
the point "C" is reduced, but a large area with smaller spin mixing is clearly present, resulting 
in b\ = 6.14 x 10 -2 . For SQA along [110], we find b\ = 6.26 x 10 -2 . This leads to an anisotropy 
A = 6%, which is still large but one order of magnitude smaller than in hep Os. This observation 
is similar to the dependence of the magnetocrystalline anisotropy energy [73] and anisotropy of 
the intrinsic anomalous Hall conductivity [T7] on the symmetry of the lattice in ferromagnetic 
crystals: the cubic W crystal exhibits a fourfold rotational axis, causing SOC to contribute to A 
in fourth order. In the uniaxial hep structure, an axis perpendicular to the c-axis is only twofold, 
and SOC enters A already in second order. Nevertheless, the comparatively large anisotropy 
of spin-relaxation in W is partly a consequence of the d-states, which yield a strong directional 
anisotropy of the Fermi surface. In contrast to this, the Fermi surface of gold consists of s-like 
states and can be regarded as almost spherical. For the Elliott-Yafet parameter in Au, we find 
a value of 6? « 3.25 x 10 -2 which is comparable in magnitude to that in W and Os, but the 
anisotropy is one order of magnitude smaller than in W. Looking at the symmetry of 6^ in 
W, we recognize that it is lower than the symmetry of the lattice. And although for a spin 
quantization axis along [001] the fourfold rotational symmetry of the lattice around this axis is 
retained by ftj^, Fig. E[d), further symmetry breaking will occur for an arbitrary direction of the 
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Figure 5: (a) Contribution to the Fermi-surface average, 6|, as function of (see text). The 
solid (dashed) arrow denotes 6? for the SQA along the c-axis (in the afr-plane). The numbers 
correspond to the respective contribution of each region to the anisotropy A, leading to a total 
anisotropy A of 59%. Note, that while the anisotropy of the EYP is clearly dominated by 
intermediate region of the major contribution to the EYP for both directions of the SQA 
comes from the region with 6^ < 0.15 (96 and 60% for the SQA along the c-axis and in the 
afr-plane, respectively). In (b) and (c), the integrated Elliott-Yafet parameter 6? is shown as 
function of the SQA direction for Os and W, respectively (different color scale). Taken from [33]. 



SQA, leaving only those point-group symmetry operations of the lattice that map the SQA to 
itself, plus the inversion symmetry k — » — k. In contrast, for the integrated value &?, Fig. [5^c), 
the full symmetry of the lattice is obviously retained. 

To conclude, we underline that the spin relaxation in metals can strongly depend on the orien- 
tation of the inject ed-electron spin axis due to a corresponding anisotropy of the Elliott-Yafet 
coefficient [33]. The anisotropy is expected to be largest in non-cubic crystals, in the presence 
of extended, nested Fermi-surface sheets that are almost degenerate, resulting in delocalized 
"spin- flip hot areas" instead of singular "spin- flip hot spots". Especially critical are cases where 
the splitting is caused primarily by the spin-orbit coupling. Since there is no theoretical limit 
on the area of the nested sheets in this scenario, the anisotropy of the EYP can be in principal 
colossal, by far exceeding the values calculated and presented here for Os, Au and W. 



4 Anisotropy of intrinsic spin Hall effect in metals 

The spin current is characterized by velocity and spin polarization. Hence, the spin current 
density Q is a tensor in R 3 ® R 3 spanned by the basis vectors ® f&. For clarity we use the 
symbols f x ,f y and f z to denote the unit vectors of spin polarization while e Xl e y and e z are 
the unit vectors of velocity. In terms of the spin Hall conductivity tensor, af- (which has three 
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indices: i denotes the direction of spin current, j the direction of applied external electric field, 
and k the direction of spin polarization of the spin current), the spin current density for a general 
direction of electric field is given by 

Q = Y t <r* j e i ®f k E j . (19) 

ijk 

While the anisotropy of the AHE manifests itself in the dependency of the magnitude of the 
conductivity vector on the magnetization direction, in the case of the SHE in paramagnets there 
is no magnetization vector M to control, only the direction of the applied electric field can be 
varied. However, the spin polarization of the induced spin current depends on the direction in 
which the spin current is measured (see Fig. [6^ a)). Hence, for a fixed electric field a given spin 
polarization s is measured only in a certain direction. Thus, in analogy to Eq. ^ we may write 

Q* = Ex<t(s), (20) 

where Q s is the spin current density for spin polarization along s and cr(s) is the SHC vector. 
(We remind that while Q s lives in the space of basis vectors e x , e y and e z , the direction of the 
spin polarization s of the spin current is spanned by basis vectors f x , f y and f z ). If the magnitude 
of the SHC vector depends on the spin polarization direction s in a material, the SHE in this 
material is said to be anisotropic. The spin Hall conductivity vector and SHC tensor, in analogy 
to the anomalous Hall effect (see introduction), are related as follows: 

^(s) = \ ^2^ijSk, (21) 

ijk 

where o\ is the l-th component of the conductivity vector, s = (s x , s y , s z ) T and is the 
Levy-Civita symbol. 

In cubic systems symmetry requires that u\- — (J z xy £ijk- Thus, the SHC may be expressed in 
terms of one material parameter, Eq. (|20|) simplifies to Q s = crLE x s, and the conductivity 



vector is <r(s) = cr z y s. Since the magnitude of the conductivity vector, <j z xyl is independent of 
s, the SHE is isotropic in cubic systems. The relationship between the direction of spin current 
and the direction of spin polarization in cubic systems is illustrated in Fig. |6^a). For the spin- 
mixing parameter and the anomalous Hall effect the dependence on the direction of the SQA 
and magnetization in the sample, respectively, can be more complicated even in cubic crystals, 
see e.g. Fig. [5^c). 

Let us consider now rigorously the situation of the SHE in transition metals with hep structure, 
Fig. [6](b), keeping in mind that the following results remain valid also for general uniaxial 
structures. If the electric field is applied along the x-direction, the magnitude of the spin 
current in y-direction will generally differ from the one in z-direction since the x-axis exhibits 
only 2- fold rotational symmetry. The spin current density in direction n = (0, cos#, sin#) T is 

n Q = -(a z xy { z cos 6 - a^Jy sin 6)E X . (22) 



Note that according to Eq. (|19j) n • Q is a vector pointing in the direction of spin polarization. 

^zy — &x y — 



We define the anisotropy of the SHE for spin polarization in the ?/z-plane as A zy = a z — a z 



Physically, if a* and a yz have the same sign, parameter A zy quantifies the difference in the 
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Figure 6: (a) Spin currents in cubic systems induced by an electric field along the z-axis. For 
electrons e with velocity v in ^/-direction the spin s points in x-direction, while for electrons 
going in x-direction the spin points in minus ^/-direction, (b) Hexagonal hep structure of the 
transition metal Ti. The spin current in direction li = (0, cos#, sin#) T induced by an electric 
field in x-direction is not perpendicular to the velocity v for a general angle i.e., v and s 
enclose an angle a ^ 90°. (c) Decomposition of the SHC of Sc into perpendicular and parallel 



components following Eq. (23). The angle a enclosed by the direction of the spin current and 
the direction of the spin polarization is also shown. At the angle #o— 62.2° the component of the 
spin polarization perpendicular to the spin current vanishes and a=180°. Taken from [29]. 



magnitude of the spin current measured along y and along z axes, when the electric field points 
along x. For a general angle 9 the components of the spin current with spin polarization parallel 
to li (Q||) and spin polarization perpendicular to n (Q±) are given by 

Qll =h'Q'h = ~A zy sm(20)E x , 

1 (23) 

Q± = (° y zx + A^cos 2 6)E X . 

If A zy 0, the spin polarization is perpendicular to li only if li is along the y or z-direction, 
otherwise spin polarization and direction of spin current enclose the angle a = arctan((5±/(5||) ^ 



90°, as shown in Fig. Mb). It follows from Eq. (23) that Q± is zero at the angle 



6 = arccos y-a% x /A zy (24) 

if a z xy and a v zx differ in sign. At this angle #o the spin polarization and the spin current are 
collinear. This is an interesting constellation, which cannot occur in cubic systems. An analogous 
situation can also occur for the anomalous Hall effect, i.e. different sign of the AHC for two 



different high-symmetry directions of the magnetization in the crystal, as we show in section 5.1 
In the latter case, there exists a direction of the magnetization in the crystal for which the Hall 
current (spin current in SHE) and the magnetization (spin-polarization in SHE) are collinear. 
Motivated by the rotational sense of the Hall current as the magnetization direction is rotated, 
we call this effect the anti- ordinary Hall effect, see section [5T| 

The case of spin current in x-direction and electric field E = (0, EcosO, Esin9) T in the yz- 
plane is simply related to the previous one by a minus sign: The components of the spin 
current with spin polarization parallel and perpendicular to the electric field E are given by 
Q|| = 7}A zy sm(29)E and Q± = —(o v zx + A zy cos 2 0)E X , respectively. At the angle 0$, Eq. (24), 
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the spin polarization and the electric field are collinear. Thus, one can achieve collinearity of spin 
polarization and electric field, or collinearity of spin polarization and direction of spin current 
if a z xy and cr v zx differ in sign. For the anti-ordinary anomalous Hall effect this means that we 
can find a direction of E such that the Hall current is collinear to the magnetization, while 
transverse Hall current is not zero. 

If the electric field is applied along the z-axis, the same magnitude of the spin current will be 
measured in all directions perpendicular to the z-axis, since the z-axis exhibits 3-fold rotational 
symmetry. The spin current in direction n = (cos#, sin#, 0) T is in this case 

n Q = (<t*J x sin 9 - a\Jy cos 9)E Z . (25) 

Symmetry requires that a zx — a yz . Consequently, the magnitude of the spin current is indepen- 
dent of 9 and the spin polarization is perpendicular to both the electric field and n. 



In the case of the hep structure the conductivity vector and the spin current density, Eq. (20), 
may be expressed in terms of the anisotropy as 

^(s) = ^s + (0,0,A^) T , ^ 
= a x yz E x s + A zy s z (E y , -E X1 0) T . 

Hence, only two parameters, a yz and A zy , suffice to describe the SHE in hep nonmagnetic metals. 
The fact that one needs only two parameters to reconstruct the exact analytical dependence of 
the spin polarization on the direction in which the spin current is measured is a manifestation 
of the geometrical anisotropy of the SHE. This is a major difference to the anomalous Hall effect 
and spin relaxation, for which the conductivity vector and the EYP have to be recalculated anew 
for each direction of the magnetization and SQA, since the EYP and the AHC for a general 
direction of the SQ A/magnetization cannot be related to the corresponding values for the high- 
symmetry axes. For example, in case of the AHE, already four parameters are needed for an 
approximate expansion of the conductivity of hep crystal up to third order in the directional 
cosines [17J. 

Next, we present in Fig. [7] the results of first principles calculations of the intrinsic SHC, Eq. ([9]), 
for the hep metals Sc, Ti, Zn, Y, Zr, Tc, Ru, Cd, La, Hf, Re and Os and for antiferromagnetic 
Cr (see section 2 for computational details). In the case of Cr we neglected the spin-density 
wave and considered the antiferromagnetic structure with two atoms in the unit cell and with the 
magnetic moments parallel and antiparallel to the z-axis. Except for Cd all metals studied in this 
work exhibit a large anisotropy of SHE, which we expect to be clearly visible in experiments. 
Of particular interest are the hep metals Sc, Ti and Ru, where the sign of the conductivity 
changes as the spin polarization is rotated from the z-axis into the xy-plane. As discussed 
before, collinearity of the spin polarization and the electric field (or the spin polarization and 
the spin current) may be achieved if the electric field (the spin current) lies in the yz-plane 
at the angle #o, Eq. (24), from the y-axis. To illustrate this we plot in Fig. [6jc) the angle 



a enclosed by the direction of the spin current and the direction of the spin polarization as 
well as the SHCs associated with Qy and Q± (see Eq. ([23])) as a function of the angle 9 for 
Sc. The critical angles at which the perpendicular component of the spin polarization vanishes 
are #o=62.2°, #o=32.1°, and 0q=19.1° for Sc, Ti, and Ru, respectively. Note that in case of 
Ru we have the case of a colossal anisotropy of the SHE: the values of the two calculated 
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conductivities differ by an order of magnitude. In the case of Cr the SHE is anisotropic as 
the cubic symmetry is broken by the staggered magnetization: If the spin polarization of the 
spin current is perpendicular to the staggered magnetization the SHC is larger by a factor 
of 1.8 compared to the case of spin polarization parallel to the staggered magnetization. We 
can thus claim that in antiferromagnets the direction of the local spins presents an additinal 
channel for the SHE anisotropy. Such anisotropy is not anymore geometrical, however, due to 
the dependence of the electronic structure on the direction of local magnetization, similarly to 
the case of the AHE in ferromagnets. 

Generally, a simple analysis of the SHC and its anisotropy in terms of a simple model becomes 
very difficult, since (i) the integrand in Eq. ^ varies very strongly as a function of k (see for 
example Fig. [9]) and the entire Brillouin zone has to be considered in the integration in order 
to reproduce the SHC quantitatively correctly; (ii) for the anisotropy of the SHC not only the 
anisotropy of the wavefunctions with respect to the SQA, discussed in the previous section, 
has to be taken into account, but also the anisotropy of the velocity matrix elements has to be 
necessarily accounted for. This makes it hardly possible to interpret the spin Hall conductivity in 
terms of a small number of virtual interband transitions. Even the sign and order of magnitude 
of the SHC are difficult to predict based on simple arguments. 

One aspect we would like to remark on is the importance of transitions in Eq. ^ which are 
driven by spin-flip SOI, and the difference between the AHE and SHE as far as the anisotropy 
of the conducitivities is concerned. Let us consider a situation of two doubly-degenerate Bloch 
states at a certain fc-point, occupied i(j n and unoccupied */j m . Let us also assume that these states 



are well-separated in energy, i.e., the first-order perturbation theory as given by equation (16), 
applies. In this case, consider the contribution to the spin Berry curvature Vl z xy which comes 
from the products of the type: 

~ (^m\Vx<J P z \^n)(^n\Vy\^m) = (<4 + + &£) {d\ + b^Vy^ + 6^), (27) 

where the SQA is chosen along the z-axis, and aln = a m \ f), etc., according to the expansion 
If we neglect the relativistic correction to the velocity operator (as our ab initio calculations show 
it is a very good approximation in most of the cases), the velocity operator does not couple states 
of different spin and the spin mixing parameter enters with the terms of the order of b m b ni which 
means that the spin-flip spin-orbit appears only in contributions to the spin Berry curvature 
which are proportional to £ 2 and higher even powers of £ (we remind that £ is the SOC strength 
in the system). It is clear that in this case the dominant contribution to the SHC comes from 
~ £ spin-conserving transitions. We prove a similar result in the next section for the anomalous 
Hall effect. 

While in the largest part of the Brillouin zone the SHC originates mainly from the spin- 
conserving SOC, in the vicinity of a degeneracy (or crossing) point, such as depicted in Fig. [8] 
(left, no SOC), both spin- flip and spin-conserving SOC can provide very large contributions 
to the SHC. Generally speaking, depending on the orbital character of the states which cross, 
the role of the spin-conserving and spin-flip SOC for the SHC around such points can be in- 
terchanged by changing the SQA. Consider for example a situation from the previous section, 
where the (doubly-degenerate) bands which cross have a dominant m and m' orbital character, 
'1^1, and are well separated from other states, Fig. [8] (left column). In this case for 
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Figure 7: For the hep metals Sc, Ti, Zn, Y, Zr, Tc, Ru, Cd, La, Hf, Re and Os and for 
antiferromagnetic Cr the spin Hall conductivities a yz and a z xy are shown as light (red) and dark 
(blue) bars, respectively. Taken from [29] . 



the SQA along z the degeneracy between the states is lifted by LS^ \ and the states keep their 
almost pure spin character. On the other hand, by pointing the SQA along x, the degeneracy 
is lifted due to LS^, and the states become strongly mixed in spin. Such anisotropy of the 
wavefunctions will contribute to the anisotropy of the SHC, but what also has to be taken into 
account is that while in the first case the v x a v z and v y velocity operators have to be considered 
in the expression above, in the second case they have to be replaced with v y a v x and v z . One has 
to realize that for a fc-point away from high symmetry directions in the Brillouin zone and for 
the bands which have mixed orbital character the mixture of LS^ and LS^ in the spin-orbit 
Hamiltonian can be complicated and the Berry curvature can vary in a non-trivial fashion with 
the SQA. Nevertheless, by looking at the situation depicted on the left of Fig. [HJ it becomes 
intuitively clear why the anisotropy of the SHE in paramagnets with structural inversion sym- 
metry is geometric: upon rotation of the SQA the spin-conserving part of SOI is continuously 
rotated into the spin- flip part, while the energy spectrum remains unchanged and the crystal 
basically remains "the same" system. Indeed, none of the gound state properties of such a 
crystal are sensitive to the direction of the SQA, and it is the non-equilibrium nature of the 
spin-relaxation and transport phenomena which makes them sensitive to it. The condition for 
such a continuous transformation between LS^ and LS^ obviously lies in the availability of 
both spin- up and spin-down states for each fc-point and energy, i.e., the spin degeneracy. 

In a ferromagnet, the spin degeneracy is lifted due to the presence of the magnetization. In this 
case, since the spin-up and spin-down subspaces are separated in energy, the spin-conserving 
SOC does not transform continuously into the spin-flip SOC at a given fc-point and energy when 
the magnetization is rotated, and two types of band crossings (or degeneracies) without SOC 
should be considered: between bands of the same, and of the opposite spin, see Fig.[8) Depending 
on the orbital character of the states, the degeneracy between them will be lifted in the first 
case by LS^ for one direction of the magnetization only (say, z), while in the second case it 
will be lifted by LS^ only for another (say, x), see Fig. pi Simply speaking, since the position 
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Figure 8: Left column: a band degeneracy due to bands of m and m! orbital character, 
\m — m'\ 7^ 1, in a paramagnet with structural inversion symmetry without SOC, is lifted due to 
spin-conserving SOC for SQA along z, and by spin-flip SOC for SQA along x. In a ferromagnet 
the band degeneracies between the bands of the same spin character (middle column), and of 
the opposite spin character (right column) have to be considered instead. Note that in this case 
the energy shifts due to SOC depend strongly on the direction of the magnetization and the 
type of the crosssing. The red (blue) color stands for spin-up (-down) character of the states, 
while green color marks the states which are of essentially mixed spin character. 

of the two types of degeneracies in energy and in the Brillouin zone, as well as their number, 
is different in a ferromagnet, effectively, for the two different directions of the magnetization we 
have two different systems with different energy spectrum and different eigenstates. In particular, 
such asymmetry is the reason for the anisotropy of the orbital moments and total energy in 
ferromagnets. Since the difference in, e.g., eigenspectrum for the two different magnetization 
directions can be hardly reconstructed analytically due to the complexity of the Halmiltonian in 
transition metals, this leads to a complicated behavior of the AHC as a function of the direction of 
the magnetization in the crystal. The same holds true for the Elliott- Yafet parameter, considered 
in the preceeding section, in which case the time-reversal symmetry in the system is effectively 
broken by a certain spin direction of the injected electron, which is able to "probe" the spin- 
mixing parameter of the states. 

Finally, we would like to make two remarks. The first one concerns the topology of the degen- 
eracy points in the Brillouin zone. The type of degeneracy shown in the upper row of Fig. [8] 
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provides the so-called "monopole" contribution to the Berry curvature [15] . intensively discussed 
in the literature, especially with respect to topological insulators [52] . Such degeneracies arise 
at single (often high-symmetry) points in the Brillouin zone, as well as along so-called "hot 
loops" \30\ YTQ . Another important contribution to the Berry curvature can be also given by 
transitions between the pairs of parallel bands degenerate along whole (often high-symmetry) 
lines or even areas in the Brillouin zone — these are the so-called "ladder transitions" [30] . In the 
language of spin-relaxation, such transitions would occur at the spin-flip hot areas in the vicinity 
of the Fermi surface. The conclusions of the discussion above hold true for both cases. Secondly, 
it is important to underline that although, referring to the perturbation theory arguments, the 
same matrix elements of SOI enter into the expressions for energy shifts, Hall conductitivities 
and spin- mixing parameter, these expressions are fundamentally different. This means that, 
for example, even though the band degeneracy in Fig. [8] would be lifted by the spin-conserving 
SOC, it can happen that the major contribution to the Berry curvature is provided by the 
spin-flip SOI, and the other way around. Concerning this issue, see also the discussion around 



second-order perturbation theory expression for the AHC in section 5.2 



5 Anisotropy of intrinsic anomalous Hall effect in metallic fer- 
romagnets 

5.1 Anisotropic AHE in uniaxial ferromagnets: first principles studies 

Experimentally, the — sometimes strong — anisotropy of the anomalous Hall effect in metals is 
a rather well-known phenomenon: see, e.g., experimental data for bcc Fe [18], fee Ni [19l [20] . 
hep Gd [21], as well as FeCr 2 S 4 [22J, Ybi 4 MnSb n [23J, YsFeiT-^Co^ [24J and R 2 Fe 17 (R = Y, 
Tb, Gd) [25J. Theoretically, the first argumentation for a strongly anisotropic behavior of the 
AHC in transition metals was provided by Roman et al [17] • In that work it was argued that 
the main reason for the observed anisotropy of the intrinsic AHC in uniaxial hep cobalt, which 
reaches as much as a factor of four between the conductivities for the magnetization in-plane and 
out-of-plane (see Table [I]), in agreement to experiment, lies in the irregular and spiky behavior 
of the Berry curvature in the reciprocal space. The main conclusions and analysis presented in 
the following three subsections holds true also for the case of the spin Hall effect. 

Let us take a look at the distribution of the Berry curvature along the high symmetry lines 
in the Brillouin zone for another uniaxial ferromagnet, Llo FePt (see Fig. [lO^c) for a sketch 
of the structure and definition of the crystallographic directions), presented in Fig. [9j The 
characteristic spikes in the vicinity of points of near-degeneracy across the Fermi energy can 
be seen, e.g., around the M-point or in the middle of the TZ-path. As we shall see in the 
next section, similarly to the case of the Elliott- Yafet parameter, the large values of the Berry 
curvature in the vicinity of such points can be inevitably related to the matrix elements of 
the spin-orbit interaction between the occupied and unoccupied states, scaled by the energy 
difference between them. As discussed for the case of the EYP and SHE, those matrix elements 
are strongly anisotropic with respect to the SQA, or, direction of M, resulting in the remarkable 
anisotropy of the Berry curvature in Fig. [9j both in magnitude (in the middle of the TZ-path) 
and sign (close to M-point). When integrated over the whole Brillouin zone, the anisotropy of 
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bcc Fe 


hep Co 


FePd 


FePt 


CoPt 


NiPt 


M||[001] 


767 


477 


135 


818 


-119 


-1165 


M||[100] 




100 










M||[110] 


810 




276 


409 


107 


-914 


M||[lll] 


842 













Table 1: Calculated from first principles anomalous Hall conductivity as a function of direction 
of magnetization in the crystal. The data are taken from Ref. [33] for bcc Fe and hep Co (similar 
values for hep Co were obtained in Ref. [32]), from Refs. [TSJ [30] for L1q FePt and FePd, and 
from Ref. [32] for Ll CoPt and NiPt. 



the Berry curvature leads to a factor-of-two reduction in the AHC in FePt as the magnetization is 
changed from out-of-plane to in-plane (see Table [TJ. In general, in cubic crystals, the anisotropy 
of the AHE with respect to the directional cosines of the magnetization appears in all odd 
(owing to the anti-symmetricity of cr with respect to M) orders starting from the third, and 
it is normally much weaker than that for the uniaxial crystals, for which the anisotropy is 
present already in the first order [17J, compare e.g. values for bcc Fe to the ones for the uniaxial 
ferromagnets in Table [T] 

When we compare the anisotropy of the AHE to the anisotropy of the Elliott-Yafet parameter 
in metals, several things come to mind. Firstly, the anisotropy of the AHC is a more complex 
quantity, which hinders analysis in terms of a simple line of arguments, as it can be done for 
the anisotropy of the EYP. This is due to the fact that while for the emergence of the EYP only 
spin-flip part of SOC plays a role and the transitions between the spin-degenerate bands can 
be ignored, for the AHC both spin-conserving and spin-flip parts of SOI have to be taken into 
account when transitions between occupied and unoccupied bands of various orbital and spin 
character varying in the Brillouin zone have to be considered according to Eq. Q for the AHC 
(see also considerations at the end of section [4]). We analyze this in more detail in the following 
two sections. Secondly, in addition to the anisotropy of the SOI matrix elements, which leads 
to the anisotropy of the wavefunctions and eigenenergies, similarly to the SHE also anisotropy 
of the velocity matrix elements matters for the total value of the AHC anisotropy. Finally, the 
Berry curvature is not confined to the Fermi surface, but has a finite spread in energy. While it 
is already clear from Eq. Q, in order to further clarify this point we refer to the distribution of 
the Berry curvature for L1q FePd alloy in Fig. [9j In this plot, the presence of wide regions in 
fc-space is evident for which the Berry curvature arises due to transitions between bands well- 
separated in energy. Since in such regions the Berry curvature also displays a very anisotropic 
behavior (see also Table [T]) it seems reasonable to ask on whether there is a certain threshold in 
energy beyond which the transitions between bands can be neglected for the anisotropy of the 
AHC. 

In order to answer this question, following Roman et al [17] , we introduce the cumulative 
anomalous Hall conductivity A(u), which accumulates all transitions in Eq. Q for which e n \^ — 
5 m k is larger than huo. In the limit of uo — >► all transitions in Eq. Q are accounted for, 
and A(cj = 0) equals the full AHC. A closer inspection of the cumulative AHC presented as a 
function of energy and magnetization direction in Fig. 12 (left) for FePt, CoPt and NiPt (the 
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Figure 9: Bandstructures and distribution of the Berry curvature along high symmetry directions 
in the Brillouin zone for the magnetization along [001] and [110] crystallographic axes in L1q 
FePt and FePd. Here tyjk] marks the (only non-vanishing) component of the Berry cruvature 
vector along the direction of the magnetization [ijk]. 
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latter two also exhibit a very large anisotropy of the AHC, see Table [T]), as well as for hep Co in 
Ref . [T7] , shows that the energetical distribution of transitions which provide the anisotropy of 
the AHC is concentrated in a narrow 1 eV window around the Fermi energy. This is in contrast 
to the distribution of A(uj) for each of the magnetization directions separately, which decays 
slowly over a much larger energy scale of several eV. This means that for the anisotropy of the 
AHC the interband transitions in the close vicinity of the Fermi level, which give very large 
contributions to the AHC, are more important than for the values of the AHC themselves. In 
principle, this is in agreement to the results for the EYP we obtained in Os, Fig. [5j there, the 
contributions from the areas where the spin-mixing parameter (Berry curvature in the case of 
the AHE) was enhanced, are dominant, while for the values of the EYP themselves the regions 
with smaller spin-mixing parameter are more important, for each of the directions of the SQA, 
see caption to Fig. [5| 

Generally, the AHC, as well as its anisotropy, displays a strong dependence on the exact position 
of band (near-) degeneracies with respect to the position of the Fermi energy. We therefore 
expect a non-trivial behavior of the AHC anisotropy for a ferromagnet with a complex electronic 
structure when the Fermi energy, or other parameters, such as the lattice constant or exchange 
splitting, are smoothly varied. An example of this phenomenon can be seen in Fig. [lO^a), in 
which the AHC and its anisotropy are plotted as a function of the "band filling" in L1q 3dPt 
alloys (see section [2] for details of the calculations). In this plot we observe that when going from 
FePt to NiPt the AHC anisotropy undergoes a change in sign and large changes in magnitude. 
With grey shaded area in Fig. flO^a) the region around the CoPt alloy is highlighted, where 
both AHC for M||[001] (cr z ) and M||[110] (cr x ) change their sign. This sign change leads to the 
occurrence of two key phenomena with respect to the anisotropic AHE. The first one — the 
colossal anisotropy of the AHE — according to calculations in Fig. 10, occurs for Fe x Coi_ x Pt 
alloy with x ^ 0.1 and for Co x Nii_ x Pt alloy with x « 0.85. For these two compounds one of the 
conductivities crosses zero, which marks the complete disappearance of the intrinsic anomalous 
Hall current Jjj for one of the magnetization directions in the crystal. This is reminiscent of 
the situation for the spin Hall effect in Ru, see Fig. [7j In terms of the longitudinal transport 
within the setup of e.g. anisotropic magnetoresistance (AMR) experiment [13], the occurrence of 
the colossal anisotropy of the diagonal conductivity would results in a metal-insulator transition 
in the crystal — in case of the colossal AHE anisotropy observed in 3dPt alloys all compounds 
remain metallic for all magnetization directions, however, and retain their complicated electronic 
structure around the Fermi energy. 



For CoPt alloy the situation, depicted in Fig. 10^d-g), is completely different. Remarkably, 
cry turns to zero at 6o = 70°, which manifests the occurrence of the anti-ordinary Hall effect 
in the crystal of CoPt, discussed already within the framework of the spin Hall effect in the 
preceeding section. At this "magic" angle, the magnitude of the anomalous Hall current J# 
is almost twice larger than it is for M||z, however, due to non-vanishing a± component of the 
AHC vector, J# is aligned along the direction of the magnetization. By analyzing Fig. 10^d-g) 
we observe that the rotational sense of the anomalous Hall current is opposite to that observed 
in the ordinary Hall effect (OHE) of free electron gas. For OHE, Lorentz forces ~ [H x v] are 
acting on electrons with velocity v in the presence of magnetic field H. The resulted ordinary 
Hall current of free electrons is always perpendicular to H irrespective of its direction, opposite 
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Figure 10: Anomalous Hall conductivity of 3dPt alloys, (a), for M along [001] (a Zl open circles) 
and [110] open squares), and anisotropy, (b) (Aa tot = a z — a Xl small open circles), with 
respect to the band filling mimicked within the virtual crystal approximation, (c) Crystal 
structure of Llo FePt alloy. Small (red) spheres stand for the Fe atoms, while large (blue) 
spheres mark the Pt ions. The primitive unit cell used in the calculations is enclosed with thicker 
lines. In the text, z stands for the [001] axis, while x stands for the [110] direction in the crystal, 
(d) Anti-ordinary Hall effect in CoPt. Red circles (blue squares) denote the a\\ (<tj_) component 
of AHC, as a function of the angle 9 of the magnetization M with [001]-axis upon rotating it 
into the [110] direction, (e)-(g) depict the relative orientation of the Hall current J#, AHC a 
and magnetization M in the situation of the anti-ordinary AHE. In (c)-(e) the magnetization 
is confined to the (llO)-plane. We note, that due to symmetry considerations, the conductivity 
vector stays together with the magnetization in the (llO)-plane. Taken from [32J. 



to the situation of the anti-ordinary anomalous Hall effect, observed in CoPt. Here, turning the 
magnetization clockwise in the (llO)-plane results in an anti-clockwise rotation of Jh , with its 
value staying rather large all the time. On the other hand, in analogy to the spin Hall effect, for 
the anti-ordinary Hall effect, it is possible to find a direction of the electric field E such that the 
Hall current is perpendicular to E and M, which are, in turn, collinear to each other. Again, 
such situation obviously cannot occur for the ordinary Hall effect of free electrons. 

In the region of 3dPt alloys in the vicinity of Llo CoPt the anisotropy of the AHE manifests 
itself in crucial ways suggesting new functionalities of the AHE-based devices. In this region, 
large changes in the magnitude of the anomalous Hall current as well as relative orientation of 
the Hall current with respect to the magnetization can be easily achieved by simple reorientation 



27 



of the sample's magnetization. While the former could be used in order to e.g. tune the relative 
magnitudes of the extrinsic and intrinsic anomalous Hall signal [TSJHQ], among most straight- 
forward applications of the latter could be a realization of the planar Hall effect (PHE) [76j [13] 
which is related to the Hall effect in ferromagnetic materials with electric field, magnetization 
and the Hall current sharing same sample plane. So far, it is believed that in most of the cases 
the PHE originates from anisotropic magnetoresistance in metallic ferromagnets, although a 
PHE mechanism stemming from the anomalous Hall effect due to non-collinearity of the mag- 
netization in semiconductor-based materials has been also suggested [77]. Within the scope of 
the anti-ordinary Hall effect, described in this work, it would be possible to observe the PHE 
coming solely from the anisotropic nature of the collinear ferromagnetic materials. 



5.2 Perturbation theory treatment: FePt 



In this section we present the perturbation treatment of the intrinsic AHC given by the linear 
response Kubo formula Eq. Q. According to first order non-degenerate perturbation theory, 
the perturbed wavefunction originating from the unperturbed wave function with spin 
a upon including the spin-orbit interaction is given by 



<j 

io = Ko) + E ■ 



£L-S 



m,0 



(28) 



where E m $ and e p ^ denote the unperturbed eigenenergies, and fc-point indices have been omitted 
for simplicity. Compared to Eq. ( fl6| ), here we consider the complete SOC Hamiltonian, and not 
only its spin- flip part. Following Eq. Q, in order to obtain an expression for a Zl the key is to 
evaluate the imaginary part of the following product: 
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(1) 



e (iK 2 



(29) 



In this expression v y and v z operators have to be considered instead of v x and v y , respectively, if 
a x is to be evaluated. The energies and e£} stand for the first-order perturbed eigenvalues. 



Substituting Eq. (28) into Eq. (29), we can sort out the terms which appear at different orders 



with respect to the SOI strength £. The purpose of this is the general analysis of simplified 
expressions, and discussion of the orders with respect to £ and their energy scales, which remain 
valid also when the (degenerate) perturbation theory is applied rigorously in higher orders. In 
the following we assume that the velocity operator does not contain the relativistic correction 
due to spin-orbit coupling, which we generally find to be a very good approximation. 



A typical first-order in £ contribution to Eq. (29) involves a sum over additional transitions via 
auxiliary states ipf^, and looks like: 



(e W -£ (1) )2 
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(30) 



while the second-order contribution to the product of the matrix elements of the velocity oper- 
ators involves already two sums of additional transitions via auxiliary states and and 
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consists of terms with the following structure: 
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m,0 
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E 



L S 



n,0 



£n,0 — £l,0 



vPmfi v y W,0 



(31) 

For the first order terms, Eq. (30), the initial state ^m0 an< ^ nna l state ^ must have the 



same spin, since the velocity operator does not act on the spin part of the wavefunction (if we 
neglect the relativistic correction to the velocity operator). This means that the state x/jf^ has 
to be of the same spin as states n and m. This can happen only due to the spin-conserving 
part of the spin-orbit interaction LS^ , as was also found by Cooper [75], meaning that only 
spin-conserving SOI contributes to the AHC in the first order with respect to £. Thus, within 
the non-degenerate perturbation theory, we would expect the largest contribution to the AHC 
from the spin-conserving part of the spin-orbit interaction. It should be kept in mind, however, 
that in materials containing heavy atoms the SOI cannot be treated as a small perturbation. 
Moreover, as follows from our previous discussion, the important role for the AHC of near 
degeneracies across the Fermi level cannot be denied, for which the above arguments, based on 
non-degenerate perturbation theory, do not apply (see also discussion at the end of section 

On the other hand, the spin-flip processes contribute only in second- and higher-order terms. We 
have also come to this conclusion for the spin Hall effect in paramagnets, following a somewhat 
different argumentation. Analysing the terms second-order in £ which come from the first-order 
perturbed wavefunctions, we find four types of contributions to the AHC, analogous to the one, 



given by expression (31). All four terms include two summations over auxiliary states ^ 



pfi 



and ^fo ' an d include two products of the matrix elements of the SOI and components of the 
velocity operator. Assuming for simplicity that the occupied state ip°' has a r — f, the four 
types of contributions can be related (omitting the energy denominators for simplicity) to the 
products of the velocity and SOI matrix elements arranged in the way, presented in Fig. [TT] In 
this figure, the diagram on the left stands for the product of ol^l^m o)> (^™ ol^^l^o)' 
(^oW^po) an d (^pol^^l^lo)' wnne the diagram on the right side stands for the product of 
<<oH<o>, ^U L S n \^i,o), (^m,oKKo> and (^ol^l^l.o)- The other two contributions 
to the second order AHC come from the diagrams in Fig. [TTJ in which the v x and v y operators 
are interchanged, while all directions of the arrows are reversed. 



As it is evident from Fig. 11, all diagrams contributing to the second-order AHC include the 
matrix elements of the spin-non-conserving part of the spin-orbit interaction. Interestingly, al- 
though a sine: le act of LSfi on a wavefunction is to flip its spin, in addition to the contribution 
to the AHC from the occupied n and unoccupied rn states of different spin character (right 
diagram), there can also be a non- vanishing contribution from the second-order transition be- 
tween the n and m states of the same spin. Finally, in addition to the contributions to the 



second-order AHC depicted in Fig. [TTJ one could in principle consider a situation when all the 
/, m-states have the same spin character and the SOI matrix elements are coming from its 
spin-conserving part. Assuming the elements of the velocity operator (at a given fc-point) to 
be real, however, it is easy to see then that the product of the SOI matrix elements would give 
a real number, since the LS^ -operator is purely complex, and thus the imaginary part of the 
second-order AHC contribution would vanish in this case. As confirmed by our first principles 
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Figure 11: Diagrams depict the second-order in SOI strength £ contributions to the AHC upon 
expanding the wavefunctions according to the first-order non-degenerate perturbation theory, 
see Eq. ( [311 ). The solid (dashed) lines stand for the matrix elements of the velocity operator 
(LS^-operator) between the bra- and ket-states, marked with indeces next to the dots, with 
the direction of the lines from the bra- to the ket-state. The horisontal thin lines separate the 
^-states from the ^-states. 



calculations, this assumption is valid. 

By calculating the spin-orbit matrix elements in the basis of states unperturbed by SOI as 
described in section 2, we applied the non-degenerate first-order perturbation theory in the 
wavefunctions with respect to SOC and computed the corresponding orders of contributions to 
the AHC in FePt. The results of these calculations are shown in Table [2] for M||x and M||z. 
Here, we considered separately the spin-conserving and spin-flip parts of our first principles 
Hamiltonian, converged the system, and applied the first order non-degenerate perturbation 
theory in wavefunctions, as described above. If we keep only LS^ spin-orbit in our calculations, 
we arrive at the value of the AHC which we denote as cr^, while keeping exclusively the LSfo 
SOC leads to the value of the conductivity a^. If and are analogously calculated non- 
perturbatively from first principles, see right column of Table [2j then, as our calculations show 

a « + a 1 ^, (32) 

where a stands for the total AHC calculated with the complete SOC Hamiltonian. Moreover, 
if the first order perturbation theory in wavefunctions is applied to evaluate the corresponding 
conductivities, the above decomposition is exact. To obtain the perturbation theory values in 
Table [2] we used a tolerance parameter A of 50 meV: that is, when for a considered unperturbed 



state \m) the difference in energy \e m ^ — s p $\ was less than A in Eq. (28), the projection on 
state \p) was considered to be zero and the corresponding term in Eq. (28) was neglected. 

Overall, by inspecting Table [2] we can conclude that the agreement of the perturbation theory 
results with first-principles results presented in Table [T] and Table [2] for FePt is reasonable. The 
major contribution to comes in the first order with respect to £ and constitutes around 
550 S/cm for both magnetization directions, while within the next contributing order these 
values are corrected by at most 30%. Thus, up to considered orders, the calculated from the 
perturbation theory agrees qualitatively with the correspoding values computed from the Berry 
curvature, and also displays a small anisotropy of the AHC when the magnetization direction 
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Table 2: Decomposition of the AHC of FePt into contributions of different orders and their 
sum (E) based on a perturbative treatment of the spin-orbit interaction, in comparison to first- 
principles non-perturbative values. All values are in S/cm. See main text for details. 



is varied, in agreement with the non-perturbative calculations. As far as the spin-flip AHC is 
concerned, the perturbation theory predicts its much smaller values as compared to the spin- 
conserving AHC. Moreover, as we can see from Table [2J the sign of is positive for M||z and 
negative for M||x, and the flip- AHC anisotropy reaches a large value of around 400 S/cm, as 
compared to the first-principles value of 320 S/cm. All of these observations are in agreement to 
the results for cr & calculated from the non-perturbative Berry curvature. A reasonable agreement 
of perturbation theory values with the first principles ones is certainly coincidencial for FePt, 
in which the AHC seems not to be dominated by transitions between the bands separated by 
less than 50 meV in energy. Nevertheless, while within our approach the singular contributions 
to the Berry curvature (arising from (near) degeneracies between the states), which can play an 
important role for the AHC, are not taken into account, certain general features of the spin- flip 
and spin-conserving AHC apparent from this analysis are universal, as discussed in the next 
section. 



5.3 Spin-flip and spin-conserving transitions 

As discussed in the previous section, and confirmed by explicit calculations within the first order 
non-degenerate perturbation theory, the spin-conserving part of SOI contributes in the first and 
higher odd orders with respect to the spin-orbit strength, while the spin-non-conserving part of 
SOI leads to non-vanishing AHC in second and higher even orders. Such oddness and evenness 
of and with respect to £ can be also demonstrated in higher orders of (degenerate) 
perturbation theory, although we do not provide explicit expressions here. As follows from 
our calculations, this remains true even when the anomalous Hall conductivity is treated non- 
perturbatively within the first principles methods, outlined in section 2. Here, we present explicit 
calcualtions of the and in L1q FePt and NiPt as a function of the spin-orbit strength in 
the system, £. The results of the calculations, in which the SOC strength was scaled uniformly 
on 3d and Pt atoms with respect to the unsealed values £0, are presented in Fig. [l3j Indeed, 
from this plot we observe that, within the accuracy of the calculations, in NiPt the conductivity 
a ^(0 — ~~ cr^(— Oj wnne m FePt cr^(£) = cr^(— £), allowing thus for an expansion of cr^(£) 
m °dd (even) powers of £. Note also that for larger values of £ the behavior of and 
a ti [ s manifestly different from linear and quadratic, respectively, marking thus the importance 
of higher order terms. 



In the latter Fig. 13, by zooming into the region of very small £, we observe that the behavior of 
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Figure 12: Total (left, black), spin-conserving (right, red) and spin-flip (right, blue) cumulative 
AHC as a function of energy in FePt, CoPt and NiPt for two different magnetization directions 
(solid and dashed lines). 




Figure 13: Left: Scaling of the spin-conserving AHC in NiPt with respect to the strength of 
spin-conserving SOC on both Ni and Pt sites. Right: Scaling of the spin-flip AHC in FePt with 
respect to the strength of the spin-flip SOC. SOC on both Fe and Pt sites was scaled uniformly. 
Inset displays the zoom-in into the region of small SOI strength. Note, that for the sake of a 
numerical experiment we also consider the regime of negative values of the SOI strength in our 
first principles calculations. 
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a tot 






p e tot 


FeTT 


Fe ri Pt tot 


PtTT 


Pt tl 


M||[001] 


818 


577 


133 


14 


18 


-27 848 


541 


282 


M||[100] 


409 


585 


-184 


210 


254 


-38 65 


426 


-361 



Table 3: AHC in FePt for two magnetization directions, resolved into spin- flip and spin- 
conserving contributions from the SOI on each atomic species (left part), compared to the 
corresponding total values (right part). All values are in S/cm. Taken from [30J. 



in NiPt (not shown) and of in FePt (inset of Fig. 13) is clearly dominated by the linear 
and quadratic in £ terms, respectively. In general, the fact that the spin-conserving transitions 
appear already in the first order with respect to £, while the spin-flip transitions in second and 
higher, has two essential consequences. Firstly, it means that the energetic spread of the spin-flip 
conductivity will be much narrower than that of the spin-conserving AHC, due to the higher 



power of the energy denominator in Eq. (131 to of the order of {s n — <s m ) , as compared to that 



of the order of {e n — e m ) in Eq. (30). This can be clearly seen in Fig. 12 (right), in which the 
cumulative AHC for FePt, CoPt and NiPt is decomposed, analogously to the total AHC, into 
spin-conserving and spin- flip contributions: 

A(u)^A^(u) + A^(u). (33) 

Noticably, while (u) decays on the scale of the bandwidth of several eV, the spin-flip cumu- 
lative AHC is localized in a much narrower energy region of the order of 1 eV. It is important to 
mention that the anisotropy of the total AHC can present a competition between the anisotropy 
of the spin-conserving and spin-flip parts, depending on the exact details of the electronic struc- 



ture, see for example Fig. [12| (right) and Fig. 10 ^b) — we refer here also to the discussion at the 
end of section HI 



Remarkably, the energetic scale of the spin- flip transitions in Fig. 12 roughly corresponds to 
the energy scale of the spin-orbit interaction of Pt atoms. This observation brings us to the 
second conclusion that we can make out from the perturbation theory analysis: the contribution 
of spin-conserving transitions to the AHC is normally dominant over the spin-flip transitions 



(this can be clearly seen in Fig. 12), since the latter appear only starting from the second order 
in SOC strength. Correspondingly, in order to promote the spin- flip contribution to the AHC, 
the spin-orbit strength in the material has to be enhanced. Let us consider this point in detail, 
and prove that the spin-flip processes in FePt are induced mostly by the strong SOI on the Pt 
atoms [30]. To do this, we selectively turn off the SOI on each atomic species inside the crystal. 
The atom-resolved spin-orbit Hamiltonian reads 

#so = £ Fe L Fe • S + £p t L pt • S, (34) 

where L M is the orbital angular momentum operator associated with atomic species n (Fe or 
Pt), and £ M is the spin-orbit coupling strength averaged over valence d-orbitals. In FePt we find 
£p e = 0.06 eV and £p t = 0.54 eV, where ^ denotes the value calculated from first-principles. 



Let us recalculate now the AHC after setting to zero either ^Fe or £pt m Eq. (34), and then 
further decompose the conductivity into the spin-flip and spin-conserving parts. The results 
are presented in Table [3| Although such a decomposition is not exact, it reproduces the total 
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Figure 14: Dependence of the total (<r* ot and <r* ot ) and spin- flip (of and of) AHC in FePt 
alloy on the strength £p t of the SOI inside the Pt atoms. Taken from [30J. 

values rather well. Namely, the sum of the total conductivities driven by SOI on Fe (Fe tot in 
Table [3]) and on Pt (Pt tot in Table [3]) is in reasonable agreement to the values of a tot for both 
magnetization directions. Moreover, the decomposition of the total atom-resolved AHCs into 
spin-conserving and spin- flip parts is almost exact, as can be seen from Table [3} Consider first 
the AHC driven by ^Fe- F° r both magnetization directions the spin- flip contribution is very 
small, while the spin-conserving part is small along [001] but large along [100]. As for the AHC 
induced by £p t , the spin-conserving part is large but fairly isotropic, while the spin- flip part is 
highly anisotropic, changing from a large positive value along [001] to a large negative value 
along [100]. This confirms that the large and strongly anisotropic o^ is governed by the SOI 
inside the Pt atoms. 

Let us confirm the conclusion we draw from the perturbation theory description via nonpertur- 
bative calculations where we tune by hand the SOI strength £p t on the Pt atoms. The results 



for the total and spin- flip AHC are shown in Fig. 14 as a function of £pt/£pf K can be seen 
that for £p t less than £p t /2, the absolute value of the spin- flip AHC does not exceed a modest 
value of 50 S/cm. In this regime a* ot and a* ot are dominated by spin-conserving processes. 
Moreover, we note that while the decrease in cr* ot is almost perfectly linear, a* ot stays fairly 
constant over a wide region of £p t values. This can be understood from the fact that for M||x 
the spin-conserving and spin-flip contributions arising from Pt largely cancel one another (see 
Table [3J, so that the total AHC is mostly driven by the SOI on the Fe atoms. In contrast, 
for M||z it is the SOI on the Pt atoms which dictates the AHC. The artificial tuning of £p t 
performed above describes rather well what happens if the Pt atoms are replaced with Pd, to 
form the isoelectronic FePd alloy [75]. This can be seen by comparing the values of a tot and 
o^ of 135 (276) and 24 (62) S/cm for M||[001] (M||[100]), respectively, in FePd with the values 



taken from the shaded area in Fig. 14, where £p t ~ £p d = 0.19 eV. In particular, the sign of the 
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AHC anisotropy in FePd, which is opposite from that in FePt, is correctly reproduced by the 
scaled calculations on FePt. 



6 Outlook 



In this review, we briefly outlined the recent progress in understanding and predicting the 
anisotropy of the spin relaxation and intrinsic anomalous and spin Hall effect in metals from first 
principles. In case of the spin-relaxation this anisotropy is the consequence of the anisotropy in 
the wavefunctions upon changing the spin-quantization axis in the crystal, which can be probed 
via a non-equilibrium process such as an injection of an electron with a certain direction of spin 
polarization into a material which exhibits the anisotropy of the Elliott- Yafet coefficient. In case 
of the spin Hall effect, in addition to the anisotropy of the wavefunctions with respect to the 
SQA, the anisotropy of the velocity matrix elements comes into play in non-cubic crystals, which 
leads to an anisotropic correlation between the direction of an applied electric field, direction 
of the spin current and its spin polarization. For ferromagnets exhibiting the anomalous Hall 
effect, in addition, eigenvalues and wavefunctions display a very non-trivial dependence on the 
direction of magnetization in the crystal, which results in a complicated relation between the 
orientation of magnetization and direction of the Hall current, as well as its magnitude. The 
anisotropy of the spin relaxation and Hall currents in perfect crystals can be so strong that it can 
reach colossal values. For spin and anomalous Hall effects, the magnitude of the Hall current can 
be even completely suppressed via a suitable choice of the direction of the electric field and/or 
magnetization. Such strong anisotropy should manifest itself clearly in an experiment, and one 
of the purposes of the current review is to stimulate further experimental studies with the aim 
of extending the functionalities of future spintronic devices. 

The phenomena considered in this work stem from the electronic structure of perfect idealized 
solids. In an experiment, especially at finite temperatures, one inevitably faces imperfections 
in the crystalline order due to impurities or disorder, phonons, magnons etc. For the Hall 
effects, disorder in the system serves as a source of additional channels for the Hall signal 
due to so-called skew- and side-jump scattering [791 (HOI ED E2]- In the perturbation theory 
picture, any sort of effects due to impurity scattering should involve the matrix elements of the 
scattering potential with the Bloch states of the perfect crystal. And while already the Bloch 
wavefunctions in the solid, as we discussed, might display an anisotropy with respect to the 
SQA, also the complicated structure of the impurity potential, especially if it is spin-polarized, 
should exhibit strong anisotropy as well given an internal spin-orbit interaction and anisotropic 
crystal field. Recently, assuming a disorder due to point-like delta-correlated defects which do 
not have any internal structure of the potential, strong anisotropy of the side-jump contribution 
to the anomalous Hall effect in ferromagnets has been demonstrated from first principles 



(c.f. Fig. 15). In this case the calculated anisotropy is a consequence of the anisotropic electronic 
structure of the perfect crystal, and the question of the anisotropy of transverse transport due 
to microscopic details of the impurity potential, which can be treated with high accuracy from 
ah initio^ still remains open and serves as a fruitful subject for future studies. 

Finally, we would like to remark, that following the same philosophy as outlined above, anisotropy 
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(b) Ni [001] <j sj (c) Ni [110] a 




Figure 15: Angle-resolved side-jump conductivity in units of S/cm as a function of direction in 
the Brillouin zone for fee Ni for two different magnetization directions. Taken from [33J. 

of the transverse transport should be also large and experimentally observable for other effects 
driven by spin-orbit interaction, such as anomalous Nernst effect in ferromagnets [83| IHlj and 
spin Nerst effect in paramagnets [85, 86J. On the other hand, in compounds which exhibit non- 
collinear magnetic order, the interplay of magnetism and spin-orbit interaction becomes very 
complex, since the non-collinearity of the local spins can effectively play the role of the spin-flip 
part of the spin-orbit interaction, and the magnetic ground state itself can be very sensitive 
to the matrix elements of the spin-orbit interaction. In such a situation, strong anisotropy of 
the transverse effects observed in this type of systems, such as magnon Hall effect [87] and 
topological Hall effect [881 ESI EQ] , is guaranteed. We are aware only of a single work in this 
direction [91j , while the phenomena mentioned above still remain largely unexplored. 
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